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Abstract In this work we present a modified TOV equation which incorpo- 
rates the cosmological constant with regard to the recent astronomical obser- 
vations that the Universe is in a phase of accelerated expansion. Using this 
modified TOV equation we considered the structure of a neutron star in SdS 
space-time and calculated maximum mass limit for neutron stars. 
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, 1 Introduction 

m 

Neutron stars can be included in the class of compact stars which are stars 
with high mass and small radius and their density will be very high. So we 
need to take into account effects from general relativity, like the curvature of 
space-time in the study of neutron stars. Baade and Zwicky were the first to 
propose the idea of neutron stars (Baade and Zwicky, 1934), pointing out that 
they have very high density and small radius. They also made a very impor- 
tant suggestion that neutron stars would be formed in supernova explosions. 
We know that there are two forces acting on a star, one of them is gravitation 
and the second one arises from the pressure. The opposing force of gravita- 
tion will be thermal pressure for ordinary stars and degeneracy pressure for 
white dwarfs and neutron stars. Thus we can obtain two coupled differential 
equations for mass and pressure which tell their change with the radius of the 
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star. Tolman-Oppenhcimcr-Volkoff (TOV) equation (Tolman, 1939; Volkoff 
and Oppenheimer, 1939) derived from the Einstein's field equation gives the 
pressure gradient of neutron stars. Numerical model of a neutron star with a 
relativistic mean field theory was given by Walecka (Walecka and Chin, 1974; 
Walecka, 1975). 

Our Universe can be represented by SdS metric, since the recent astro- 
nomical observations show that the expansion of our universe is accelerating 
(Guzzo et al., 2008; Perlmutter et al., 1999; Riess et al., 1998; Weinberg, 
1972). Thus Schwarzschild-de Sitter metric described with a positive cosmo- 
logical constant and having two horizons is considered to be a model of this 
universe. In such a case, the use of TOV equation, derived from the Einstein's 
field equation in the case of accelerating universe, will be more apt for the 
theoretical constructions of stellar models. Thus we will get a modified TOV 
equation for pressure, for calculating the maximum mass of a neutron star. 
Anisotropic models for compact self gravitating objects have been studied ex- 
tensively (Corchero, 2001; Dev and Gleiser, 2003; Hcrrera and Santos, 1997; 
Ivanov, 2002; Mak and Harko, 2003). Our aim is to establish a modified TOV 
equation which is the prime structure equation for compact stars. 

In §2 we will show the derivation of the new equation to replace the TOV 
equation, starting with the standard metric for a spherically symmetric star 
and the consecutive use of Einstein's field equation in the case of accelerating 
universe. Then we will show the numerical calculation of maximum mass of 
Neutron stars by plotting the mass-radius relation. 



2 Derivation of the modified TOV equation 

We have the standard metric (Weinberg, 1972) for a spherically symmetric 
star, 

ds 2 = -U(r)dt 2 + V(r)dr 2 + r 2 d6 2 + r 2 sin 2 6dcf , (1) 

where, 
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FW " WY (2) 



Ricci tensor is given by, 



i,h dx k 
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Substituting for the Christoffcl symbols, we will get the equations as, 



"~ 2U AU\V + U) rV 
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V + U ) + V 
R(jxj> ~ R00, (6) 



R " = - 1 + W (-v + u + y' (5) 
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U" U' (V U'\ U' 
^ = -777 + 777 77 + 77 ~ 777" ( 7 ) 



2V W V V U J rV 

The result ~ Rgg is the consequence of the rotational invariance of the 
metric. 

Einstein's field equation in the case of accelerating universe (Weinberg, 
1972) is given by, 

G^v = -&rGTp, + Ag„ v , (8) 

where 

G>„ = R^u — -g^R. (9) 
Substituting Eq. © in Eq. ©, we get 

R MV = -SttG (t^ - \g^Tl - ^g^ ■ (10) 

From hydrodynamics, we have the relation for energy-momentum tensor in 
the case of models with cosmological constant, 

Tftv = pg^v + (p + e) UpUv, (11) 

where p-the modified pressure, e-the modified total energy density and u^-the 
velocity four vector. Here 

P = p -^G i = e + ^G> (12) 
Also the velocity vector is defined so that, 

g^u^Uu = -1. (13) 

Assuming fluid at rest, we have 

u r = tig = it0 = 0, (14) 

u t = -(-g a r h = -Vm- as) 

Our assumptions of time independence and spherical symmetry imply that 
p and e are functions only of the radical co-ordinate r. 
Using Eq. (TT0|) , we will get 

R -E1-EL(YL + EL\ -IYL 



2U AU V V U J r V 

= -4nG(e-p)V-AV, (16) 
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= -47rG(e - p)r 2 - Ar 2 , (17) 
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-AttG(€ + 3p)U + AU, 



(18) 



Here single prime denotes first derivative and double prime denotes second 
derivative. 

We have the equation for hydrostatic equilibrium, 



U' 



-2p> 



U p + e 



(19) 



Let us derive the equation for V(r) alone. Using Eq. (fH)|) , Eq. (fT7|) and Eq. 
(fig)), we will get 

-V 1 1 



InGe - A, 



(20) 



which can be written as 



= 1 - 8nGer 2 - Ar 2 . 



Integrating, 



V(r) 



1 - 



2GM(r) A 2 
r 3 



where 



M(r) = / 4irr' 2 e(r')dr' 



(21) 



(22) 



(23) 



Eq. (|2"2"|) tells that outside the star, the space-time looks like SdS. 

Now, we can use Eq. (fl9"]h Eq. (J22J) and Eq. (J23J) to eliminate U(r) and 
V(r) from Eq. p?|) . which becomes 



e(r) 



r 2 p' = GM (r) e (r) 



Thus the modified TOV equation is, 



dp 
dr 



Arrr p(r) , ylr 3 
M(r) GM(r) 



2GM(r) _ Arf_ 
r 3 



(24) 



GM (r) e (r) 
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(25) 



This equation is different from TOV equation and this becomes TOV equa- 
tion for A = 0. 

Now we can use this new equation to find the maximum mass of neutron 
stars. In this work, neutron stars are considered to be composed of a Fermi 
gas of degenerate neutrons. The mass-radius and density-radius relations are 
plotted and maximum mass limit is checked. 
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We have the standard structure equations for a spherically symmetric star 
which describe how the mass and the pressure of a star change with radius: 

dm _ Anr 2 e{r) 
dr c 2 ' 

dp _ Ge{r)m(r) 

t 2 2 ' \ ') 

dr c^r^ 

where G is Newton's Gravitational constant, e is the mass density, m is 
the mass up to radius r and c is the velocity of light. 

The modified TOV equation is similar to the differential equation for the 
pressure, Eq. (|27|) . but here has three correction factors. All three correction 
factors are larger than one, i.e. they strengthen the term from Newtonian 
gravity. Here one has to solve (Chandrasekhar, 1931) the coupled differential 
Equations (|26|) and (|25jl . We must notice that while ^ is positive, $ must 
always be negative. Starting with certain positive values for m and p in a small 
central region of the star the mass will increase while the pressure decreases 
eventually reaching zero. We can set m(r = 0) = and , in addition we specify 
the central pressure p{r = 0) = po in order to solve for Equations (|26[) and 
(I25p . The behaviour of m and p as a function of radius will become important 
for numerical calculations. 

Both the coupled differential equations depend on energy density. So we 
need a relation between pressure p and energy density e which is known as the 
Equation of State(EoS) (Shapiro and Teukolsky, 1983) for the neutron stars. 
For mass calculation, we will use the well known EoS: 

p = Ke\ (28) 

where K is a constant and 7 is the polytropic exponent. 7=3 in general rela- 
tivists case (Silbar and Reddy, 2004). 

We will give the value of cosmological constant (Carmeli and Kuzmenko, 
2001; Carroll, 2000) as A^10~ i6 kmr 2 . 



3 Result and Discussion 

The mass-radius and density-radius relations obtained is shown in Fig. [T] and 
Fig. [2] For calculational purpose, we put G=l and c=l and then pick the units 
of all dimensionful quantities to be powers of kilometers. 

Converting the results back from natural units, we will get the maximum 
mass of neutron star M=0.74M Q and radius R=6.98km. 

Comparing with the result that obtained when calculated using the nor- 
mal TOV equation, we found that the obtained values are lesser. So we can 
conclude that there is a slight influence of cosmological constant in the stellar 
structure. The influence of cosmological constant will be higher if the value is 
higher than the recently accepted one. Future advancements in dark energy 
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Fig. 1 Mass-Radius plot for neutron stars for a central density 1 X W le gm/cm 3 
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Fig. 2 Density-Radius plot for neutron stars for a central density 1 X 10 lfi gm/cm 3 



research may strengthen the importance of cosmological constant. The modi- 
fied TOV equation could also be used in the calculations of surface tension of 
compact stars (Bagchi, et al., 2005; Dey, et al., 1998). 
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